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Unstructured grid large eddy simulation 
of wall bounded turbulent flows 

By Kenneth Jansen 
1. Motivation and objectives 

Historically, large eddy simulations (LES) have been restricted to simple geome- 
tries where spectral or finite difference methods have dominated due to their efficient 
use of structured grids. Structured grids, however, not only have difficulty repre- 
senting complex domains and adapting to complicated flow features, but also are 
rather inefficient for simulating flows at high Reynolds numbers. The lack of ef- 
ficiency stems from the need to resolve the viscous sublayer which requires very 
fine resolution in all three directions near the wall. Structured grids make use of 
a stretching to reduce the normal grid spacing but must carry the fine resolution 
m the streamwise and spanwise directions throughout the domain. The unneces- 
sarily fine grid for much of the domain leads to disturbingly high grid estimates. 
Chapman (1979), and later Moin & Jimenez (1993), pointed out that, in order to 
advance the technology to airfoils at flight Reynolds numbers, structured grids must 
be abandoned in lseu of what are known as nested or unstructured grids. Figure 1 
illustrates the ability of an unstructured mesh to refine only the near wall region. 
Note the large number of points near the wall (where the fine vortical features need 
better resolution) and the coarseness in all directions away from the wall (where the 
scales are much larger). The important difference between this approach and the 
usual structured grid stretching is that the number of elements used to discretize 
the spanwise and streamwise features of the flow is reduced in each successive layer 
coming off the wall. This is due to the fact that the elements not only grow in the 
norinal direction, but in the other directions as well. This greatly reduces the total 
number of points or elements required for a given Reynolds number flow. 
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Figure 1 . An unstructured grid places a large number of points at the wall but 
remains coarse in the freestream. The full mesh is shown in part (a) while a zoom 
of near wall corner is shown in part (b) to illustrate the refinement. 
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To estimate the number of unstructured grid elements required to simulate an 
airfoil, Chapman used flat plate skin friction analogies and a computational domain 
extending one-fifth chord in the spanwise direction to obtain 


N = 0.2 


Re^ 

A+A+ 


( 1 ) 


where A+ and A+ are the grid resolutions in the streamwise and spanwise directions 
on the body surface and Re c is the Reynolds number based on the chord, c and the 
freestream velocity, u in f. This estimate assumes that the fine resolution near the 
wall is carried out for 10 layers to accurately resolve the viscous sublayer. Then, 
outside of the viscous sublayer, the elements grow rapidly in all three dimensions 
with increasing distance from the wall as described above. Moin k Jimenez suggest 
that current subgrid-scale models should allow A + = 200 and A+ = 50. When these 
values axe substituted into (1), we observe that approximately 1.2 x 10 elements 
will be required for airfoils with a chord Reynolds number of Re c = 10 and 80 x 10 
elements for the more practical flight Reynolds number of Re c = 10 . Simulations 
of this scale are possible on today’s supercomputers. 

The use of unstructured grids, coupled with the advances in dynamic subgrid- 
scale modeling such as those made by Germano et al. (1991) and Ghosal et al. 
(1992), make LES of an airfoil tractable. The finite element method can efficiently 
solve the Navier-Stokes equations on unstructured grids. Although the CPU cost 
per time step per element is somewhat higher than structured grid methods, this 
effect is more than offset by the reduction in the number of elements. 


2. Accomplishments 

2.1 Computer code 

The proposed finite element formulation is based on the work of Jansen et al. 
(1993), who used the method to model the compressible Reynolds-averaged Navier- 
Stokes' equations. These simulations were performed by time marching a transient 
simulation to a steady solution. The code was optimized for rapid convergence 
without regard for time accuracy. For the current work, greater attention has been 
given to the efficient time accuracy before application to LES. Both explicit and 
implicit time integration methods have been developed and tested. Currently, the 
formulation has two implicit methods (first-order for acceleration towards a steady 
state and third-order for time accurate integration) and a higher order accurate 
family of explicit time integration methods. 


2.2 Time step estimates 

There is, in general, a tradeoff between explicit methods, which are cheaper per 
time step, and implicit methods, which require fewer time steps due to the avoidance 
of the stability limits. It can be shown through methods similar to Chapman s 
spatial estimate that the viscous stability limit leads to the following time step 
limit . + 

A" - T 

1 C fm „Re c 


( 2 ) 
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where T is the time it takes the mean flow to cross the chord of the airfoil (T = 
c / u *nf ) an d Cf mmx is the maximum of the coefficient of friction. 

There is also a stability limit associated with advection. The time step associated 
with this stability limit is a little more difficult to estimate since it depends on both 
the mean flow advection, u, and the length of the element in the flow direction, A,. 

hese quantities vary throughout the flow. Assuming a logarithmic velocity profile 
and geometric stretching of the elements coming off the wall, it can be shown that 
the critical point occurs in the buffer layer near y + = 10. Respecting this advective 
stability limit leads to the following advective time step limit 


where 


A“ = aA" 

(3) 

0.6A+ 

a * A +2 

(4) 


For the problem proposed above, A+ = 200 and A+ is expected to be near 1.0, 
which makes a 120. This time step corresponds to a A+ w 120, which will not 
yield sufficient accuracy. Therefore, this stability limit is not likely to have any 
bearing on the size of the time step. 

Since we are solving the compressible Navier-Stokes equations, a third time step 

restriction must be considered. The acoustic stability limit can be estimated as 
follows 


a? = m? 


where 

8 _ 

At 


(5) 

( 6 ) 


where a is the acoustic speed and M = u M /a is the freestream Mach number. 
Clearly, (3 is less than one, making this the most restrictive stability limit. 

If A y is equal to one, it can be shown that A,' corresponds to a A = 1.0. Cur- 
rent channel flow LES simulations have had success with Af = 10.0. Assuming that 
this temporal resolution is adequate for the airfoil, the acoustic time step stability 
limit will be far too restrictive. The implication of this result is that compressible 
formulations must provide an implicit treatment of the advective term. Further- 
more, special care must be taken to show that the method is not adversely affected 
by simulating the flow at very high acoustic CFL numbers. For the conditions 
described above, typical flows lead to the following acoustic CFL estimate 


aA t 10.0 70.7 

~Al = ~*~M~ (7) 


which exceeds 350 for a Mach number of 0.2. 
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FIGURE 2. Unstructured grid for flow around a cylinder at Red — 100, M — 0.2. 
This mesh contains 4004 elements which corresponds to 2056 nodes for a linear 
space. 


2.3 Preliminary simulations 

To verify that the modifications to the code achieve time accuracy, the method 
was applied to laminar flow over a cylinder at Red = 100, M = 0.2. This flow 
leads to periodic vortex shedding and, therefore, gives some measure of a method’s 
temporal- accuracy. The unstructured triangular mesh is shown in Figure 2. Note 
the local refinement near the cylinder and in the wake. The lift coefficient obtained 
by using piecewise linear shape functions in space and time can be seen in Figure 
3. The Strouhal number for this discretization is 0.167. The acoustic CFL for this 
problem is 20.0. 

3, Future plans 

3.1 Time integration 

The cylinder problem is not an adequate test of the formulation s ability to run 
at very high acoustic CFL numbers. A channel flow at a higher Reynolds will be 
run to determine the upper limit for the acoustic CFL number. Should problems 
arise, a change to an incompressible formulation may be appropriate. Such finite 
element formulations are currently being used by Hauke & Hughes (1993) and Simo 
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Figure 3. Periodic vortex shedding illustrated through the lift coefficient. 

& Armero (1993) for laminar flows. The drawback to an incompressible formulation 
is that many interesting airfoil problems require the consideration of compressibility 
effects. 


3.2 Mesh generation 

The mesh requirements of the airfoil problem exceed existing mesh generation 
capability. In order to stretch the elements to the level described above and to reduce 
the number of elements with distance from the wall, new mesh generation techniques 
must be developed. It is crucial that elements do not have angles which approach 
180°. An algorithm to accomplish this goal for airfoils is under development in 
collaboration with Tim Barth of NASA Ames. This algorithm also should provide 
smoother element shape changes than those observed in Figure 1, resulting in higher 
quality solutions. 


3.3 Subgrid-scale modeling 

The dynamic models developed at CTR need to be implemented into the un- 
structured grid code. This is not expected to be too difficult, especially for higher 
order elements which have a built-in test grid (the corner nodes). 

3.4 Further speeding up of the code 

The code has been largely optimized for marching to steady state solutions. It 
may be possible to further optimize the code for time accurate calculations. The 
code currently runs at 440 MFLOPS on the Cray 090 and 25 MFLOPS per processor 
on the Thinking Machines CMS. These execution rates are quite fast, but some 
savings in the number of FLOPS per time step per element may be attainable. 
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S.5 Airfoil simulation 

Upon the successful completion of these tasks, the code will be applied to air- 
foil problems with the ultimate target being airfoils at or near maximum lift; see 
Coles and Wadcock (1979). These flows commonly have separation bubbles that 
are difficult to predict with Reynolds-averaged Navier-Stokes equation models and, 
therefore, present an opportunity to demonstrate the utility of LES approaches. 
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